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Abstract. We have analyzed survey and pointed ROSAT 
PSPC observations of the Seyfert galaxy NGC 3227. Large 
amplitude X-ray variability is detected, with a factor 15 
change in count rate within about 3 years. Smaller changes 
are seen on the timescale of days, the largest being a factor 
of 3.5. No strong spectral variability is found throughout 
the pointed observation. The X-ray spectrum is modeled 
in terms of warm absorption and both, a dust-free warm 
absorber and one with internal dust, give an excellent 
spectral fit. Additionally using multi- wavelength observa- 
tions of NGC 3227, the dusty warm absorber is favored 
and suggested as an explanation for the lack of a large 
cold absorbing column in the soft X-ray region despite 
excess reddening along the line of sight. Self-consistent re- 
calculation of the models including Galactic-ISM-like dust 
results in an ionization parameter log U ~ —0.25 and a 
column density log A^ — 21.8 of the warm absorber. X- 
ray variability and dust survival arguments are used to 
constrain the density and further properties of the ionized 
material. The influence of dust and other parameters on 
the thermal stability of the warm absorber is investigated. 
Predictions are given for the absorber-intrinsic optical-UV 
line emission and absorption. 

Key words: Galaxies: active - individual: NGC 3227 - 
emission lines - Seyfert - X-rays: galaxies 
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1. Introduction 

NGC 3227 is a Seyfert 1.5 galaxy at a redshift of z=0.003. 
It shows signs of interaction with its dwarf elliptical com- 
panion NGC 3226. Although studied in many spectral re- 
gions, most attention has focussed on the optical wave- 
length range: NGC 3227 has been the target of several 
BLR mapping campaigns (Rosenblatt et al. 1992, Sala- 
manca et al. 1994, Winge et al. 1995) and high-resolution 
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line-profile studies (e.g. Whittle 1985, Rosenblatt et al. 
1994). High spatial resolution observations of the kine- 
matics and ionization structure of the circum-nuclear re- 
gion were presented by Mediavilla & Arribas (1993) and 
Arribas & Mediavilla (1994). Rubin & Ford (1968) and 
Mundell et al. (1995a) reported evidence for very high 
emission-line reddening. 

X-ray variability of NGC 3227 was detected with 
HEAO 1 by Tennant & Mushotzky (1983), who found 
the flux to increase by a factor of 1.4 within 5 hours. Vari- 
ability was also found in the EXOSAT data (Turner & 
Pounds 1989) and in an ASCA observation (Ptak et al. 
1994). 

Reichert et al. (1985), using Einstein data, favored a 
partial covering model to describe the (0.8-4.5) keV X-ray 
spectrum of NGC 3227. Differential variability between 
the EXOSAT LE (0.05-2 keV) and ME (2-10 keV) in- 
struments was interpreted in terms of either partial cover- 
ing with a strongly varying cold column or a variable soft 
excess (Turner & Pounds 1989). An unusually flat pow- 
erlaw of photon index Fx ~ -1.5 was obtained. An iron 
K-emission line near 6.5 keV was detected in the Ginga 
spectrum of NGC 3227 (Pounds et al. 1989). This ob- 
servation was modeled by George et al. (1990) in terms 
of reflection from a cold accretion disk. The underlying 
X-ray continuum emission was then best described by a 
powerlaw with index = -1.86±0.04. Given the lack of 
spectral resolution below ~ 1 keV in these studies, it was 
not possible to distinguish between alternative descrip- 
tions of the soft spectral shape and to evaluate the cause 
of variability. 

Evidence for the presence of a warm absorber in NGC 
3227 was indicated by Netzer et al. (1994) by the use of X- 
ray color diagrams, and reported by Ptak et al. (1994) in 
an ASCA observation. As the cause for variability detected 
in the ASCA data, Ptak et al. found two possibilities. One 
is a change in powerlaw index by AFx ~ 0.2 and constant 
properties of the ionized material. The other is complex 
behavior of the warm absorber with variability of both, 
ionization parameter and warm column density. 
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For the present study, we have analyzed the survey 
and archival pointed ROSAT observations of NGC 3227 
performed with the PSPC (Triimper 1983; Pfeffermann 
et al. 1987). Since there is non- X-ray evidence for a large 
cold column of absorbing material along the line of sight 
(e.g. Cohen 1983), and as the absorption structure of the 
comparatively lowly ionized warm absorber extends well 
below 0.5 keV, ROSAT is particularly well suited for an 
investigation of these features. Furthermore, in a sequence 
of studies to assess the possibility of a warm-absorber con- 
tribution to one of the known high-ionization emission-line 
regions in individual AGN (Komossa & Fink, e.g. 1997b), 
NGC 3227 marks the low-ionization end of the observed 
warm absorbers. In this context, it is also interesting to 
note that Salamanca et al. (1994) reported the existence 
of broad wings in Ha in NGC 3227 that do not follow 
the continuum variability, indicating a higher than usual 
degree of ionization, as would be expected for a warm ab- 
sorber. 

The paper is organized as follows: In Sect. 2 we present 
the observations. The data are analyzed with respect to 
their spectral and temporal properties in Sects. 3 and 4, 
respectively. In Sect. 5, evidence for the presence of dust 
mixed with the warm absorber is presented and compared 
to alternative scenarios. Further properties of the ion- 
ized material are derived and discussed, and the absorber- 
intrinsic line emission and absorption in the optical-UV 
spectral region is studied. A summary and the conclusions 
are provided in Sect. 6. 

A distance of 18 Mpc is adopted for NGC 3227 as- 
suming a Hubble constant of Ho = 50 km/s/Mpc and the 
galaxy to follow the Hubble flow. 

If not stated otherwise, cgs units are used throughout. 



2. Data reduction 

2.1. Pointed data 

The observation was performed with the ROSAT PSPC 
from May 8 - 19, 1993, centered on NGC 3227. The to- 
tal exposure time is 19.6 ksec. The source photons were 
extracted within a circle chosen to be large enough to 
ensure that all of the source counts were included, but 
small enough to contain negligible contribution from the 
much weaker neighboring source NGC 3226. The back- 
ground was determined after removing all detected sources 
within the inner 19' of the field of view. The data were 
corrected for vignetting and dead-time using the EXSAS 
software package (Zimmermann et al. 1994). The mean 
source count rate is about 0.5 cts/s. For the spectral anal- 
ysis source photons in the amplitude channels 11-240 were 
binned according to a constant signal/noise ratio of 21a. 
For the temporal analysis the minimal bin size in time was 
400 s to account for the wobble mode of the observation. 



2.2. Survey data 

NGC 3227 was observed during the ROSAT all-sky sur- 
vey (RASS) from Nov. 12 - 15, 1990, with an effective 
exposure time of 610 s. It is clearly weaker in this ob- 
servation as compared to the pointing. However, due to 
the larger effective pointspread function during the sur- 
vey observations, it was not possible to exclude the con- 
tribution of the weak neighboring source NGC 3226 (al- 
though the position of the centroid of the X-ray emis- 
sion agrees with the coordinates of NGC 3227). There- 
fore, 'source'photons within a circle including both objects 
were extracted. The background was determined from two 
source free sky fields along the scanning direction of the 
telescope. The data were corrected for vignetting. A to- 
tal number of 52 source photons was detected. Assuming 
NGC 3226 to have remained constant in count rate from 
the RASS to the pointed observation (with 0.018 cts/s), 
results in a coimt rate of 0.067 cts/s for NGC 3227, re- 
vealing strong variability. Unfortunately, the low number 
of photons accumulated during the RASS does not allow 
a more detailed discussion of the spectral and temporal 
properties. 

3. Spectral analysis 

3.1. Standard spectral models 

A single powerlaw (with cold absorption column as a free 
parameter) is a bad description of the soft X-ray spectrum 
(x^red=l-7; Table 1). The resulting powerlaw is very Hat, 
with an index Fx =-1.2, and strong systematic residuals 
remain (Fig. 1). None of several single-component models 
that were compared with the data provides an acceptable 
fit. This also holds for 2-component (powerlaw plus soft 
excess) models in which the powerlaw index Fx is fixed to 
-1.9 with all other parameters free. 

A powerlaw of free index with a soft excess, or a warm- 
absorbed powerlaw of canonical index (next section) give 
successful fits. The soft excess was parameterized as black 
body, for which we find fcTbb — 0.07 keV, or the standard 
accretion disk model after Shakura & Simyaev (1973). The 
latter yields an unusually low black hole mass of Mbh — 3 
X 10^ M0 for fixed accretion rate M/Medd = 1) and a cold 
column density of ~ 0.70 x 10^^ cm'^ (Table 1). Since 
also the underlying powerlaw is again unrealistically fiat, 
with Fx — -1.1, in contradiction to the higher-energy ob- 
servations (e.g. George et al. 1990, Fx — -1.9), this model 
is not discussed further. 

Finally, we applied a partial covering model to the 
data, in which part of an intrinsic powerlaw spectrum is 
directly seen, and part is absorbed by a very high cold 
column (A^H left free or fixed to 8.3 xlO^^ cm~^; cf. Table 
2). This approach is motivated by the evidence of large 
amounts of cold matter along the line of sight deduced 
from emission-line reddening, as is further discussed be- 
low. No successful description of the data was achieved. 
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3.2. Warm absorber models 



3.2.1. Model properties and assumptions 

To model the X-ray spectrum in terms of warm absorp- 
tion, we have calculated a sequence of photoionization 
models for the warm material, using the code Cloudy 
(Ferland 1993). Wc assume the material (a) to be pho- 
toionized by continuum emission of the central pointlike 
nucleus, (b) to be one-component and of constant den- 
sity, and (c) to have solar abundances (Grevesse & An- 
ders 1989). In those models that include dust mixed with 
the warm gas, dust properties like those of the Galactic 
diffuse interstellar medium are adopted (Mathis, Rumpl 
& Nordsieck 1977; MRN) if not mentioned otherwise, and 
the gas-phase abundances are depleted correspondingly (a 
mean of Cowie & Songaila 1986), as in Ferland (1993). The 
MRN dust consists of two grain species, graphite and as- 
tronomical silicate, and a powerlaw grain size distribution. 

The ionization state of the warm absorber can be char- 
acterized by the hydrogen column density of the ion- 
ized material and the ionization parameter U , defined as 



U = Q/iAnr^nnc) 



(1) 



where Q is the number rate of incident photons above 
the Lyman limit, r is the distance between nucleus and 
warm absorber, nn is the hydrogen density (fixed to 10^'^ 
cm~^ unless noted otherwise) and c the speed of light. 
Both quantities, iV„ and U, are determined from the X- 
ray spectral fits. 

The spectral energy distribution (SED) chosen for the 
modeling corresponds to a 'mean Seyfert' contimmm con- 
sisting of (i) piecewise powerlaws from the radio to the 
gamma-ray spectral region, or (ii) an additional EUV 
bump parameterized as a black body of T = 120000 K 
contributing the same amount to Q as the powerlaw com- 
ponent. As we have verified in a series of test calcula- 
tions, the EUV spectral shape generally has no strong 
influence on the X-ray properties of the warm gas, its 
ionization state is dominated by the X-ray part of the 
incident continuum (Komossa & Fink 1997a). This also 
holds for the case of NGC 3227. (The re-calculation of 
the models to include an additional EUV bump was per- 
formed to check the robustness of predictions made for 
the absorber-intrinsic emission and absorption in lines in 
the optical and UV spectral region; cf. Sects 5.3.1, 5.3.2.) 
Furthermore, the observed UV continuum of NGC 3227 is 
unusually steep which seems to be due to strong reddening 
(Sect. 5.1.2), whereas the warm absorber most probably 
sees the unreddened continuum. More explicitly, the mean 
Seyfert SED employed consists of an UV-EUV powerlaw 
of energy index Q!uv-x=-1-4 (Kinney et al. 1991), extend- 
ing up to 0.1 keV, a mean continuum after Padovani & 
Rafanelli (1988) from the radio to the optical region with a 
break at 10/im and an energy index a = -2.5 A-longwards, 
and an X-ray powerlaw that breaks into the gamma-ray 
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Fig. 1. The upper panel shows the observed X-ray spectrum of 
NGC 3227 (crosses) and the best-fit dust-free warm absorber 
model (solid line). The second panel displays the fit residu- 
als for this model, and the third panel the residuals from the 
best-fit dusty warm absorber. For comparison, the residuals 
resulting from a single powerlaw description of the data are 
shown in the lowest panel (note the different scale in the ordi- 
nate). 
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Table 1. X-ray spectral fits to NGC 3227 (pi = poworlaw, disk = accretion disk model as described in the text, wa = warm 
absorber). The errors are quoted at the 90% confidence level. 



model 


[10^1 cm-2] 


log U 


log _/Vw 

[cm-2] 


log Normpi 


1 X 


[10* Mo] 


iVi /iViedd 


Xredl"-0-/j 


pi 


0.38±0.02 






^•"'-'-O.Ofj 


-1.19±0.11 






1.67(18) 


pl+disk 


0.70±0.16 






^•"'-'-0.12 


-1.1±0.6 


3.0±0.8 


1.0(3) 


0.58(16) 


wa 


0.55±0.13 


-1.03±0.10 


21.49±0.13 


-4.24±0.02(2) 


-1.9(3) 






0.78(17) 


dusty wa 


0.55±0.04 


-0.25±0.25 


21.75±0.10 


-4.22±0.03(^^ 


-1.9(3) 






0.69(17) 



(^) Normalization at 1 keV at 10 keV (3) fixed 



region at 100 keV. Since it turned out that the data do 
not allow to well constrain the underlying X-ray powcaiaw 
index, it was usually fixed to (i) the canonical value of Fx 
= -1.9 (used if not stated otherwise), that was found for 
NGC 3227 by George et al. (1990) in the Ginga observa- 
tion, and (ii) Fx = -1.6, the mean value derived by Ptak 
et al. (1994) using ASCA data. 

3.2.2. Model results 

Fitting the warm absorber model to the X-ray spec- 
trum of NGC 3227, we find an ionization parameter of 
\ogU ~ — 1.0 and a column density of the ionized mate- 
rial of logA/w 21.5, with xVdof = 13.3/17 (Table 1). 
The mean observed X-ray flux for this model is / = 2.15 
xlO~^^ erg/cm^/s, corresponding to an intrinsic (0.1-2.4 
keV) luminosity corrected for cold and warm absorption 
of ix = 0.8 X lO''^ erg/s. The contribution to the ROSAT 
spectrum of emission and reflection from the warm mate- 
rial, calculated with the code Cloudy for a covering factor 
of the ionized absorber of 0.5, is found to be negligible. 
The cold column, Nh ~ 0.55 x 10^^ cm~^ is larger than 
the Galactic value (iVg"' ~ 0.22 x lO^^ cm'^; derived upon 
interpolation of the Dickey & Lockman (1990) data), but 
not as large as implied by emission-line reddening, as is 
further discussed in Sect. 5.1. Dust mixed with the warm 
gas could supply the reddening without implying a corre- 
sponding cold column. 

Since the presence of dust clearly modifies the result- 
ing X-ray absorption structure (Komossa & Fink 1997a, 
1997c), we re-calculated the models, now including dust 
with Galactic ISM properties. (A density of the warm ma- 
terial of log riH = 7 was used in these models. Whereas 
usually the X-ray absorption does not strongly depend on 
density, this value was chosen to ensure dust survival.) 
Applying the dusty absorber model to the data provides 
an excellent fit with xVdof = 11.7/17 (Fig. 1). We find 
logU c± —0.25 and logA^w — 21.8 and a cold column of 
iVn — 0.55 X 10^^ cm~^ for Fx = -1.9. (For comparison, 
the corresponding values for the Fx = -1.6 model are: 
logC/ ~ -0.6, logTVw 21.5, A^H =i 0.47 x 10^^ cm'^, and 
xVdof = 13.9/17.) 



The residuals from the best spectral fits, as compared 
to the single powerlaw model, are displayed in Fig. 1. The 
unfolded X-ray spectrum is shown in Fig. 5. No distinc- 
tion between the dusty and dust-free absorber model is 
possible on the basis of the quality of the X-ray spectral 
fits (Table 1). 

4. Temporal analysis 

4.1. Flux variability 
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Fig. 2. X-ray lightcurve of NGC 3227 during the pointed obser- 
vation, binned to time intervals of 400 s. The time is measured 
in seconds from the start of the observation. 



The X-ray lightcurve for the pointed observation is 
shown in Fig. 2. With respect to the minimal resolved time 
interval of 400 s, the observed count rate varies between 
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0.27 cts/s and 0.97 cts/s, corresponding to an amplitude 
of variability of a factor of 3.5 (occuring within the to- 
tal observing period of about 10.8 days). The shortest re- 
solved doubling timescale is 380 minutes and the strongest 
change within 400 s is a factor of 1.5 drop in count rate. 

Even stronger variability is detected between survey 
and pointed observation (separated by ~ 3 y) , with a fac- 
tor ~ 8 variability in the mean count rates, and a factor of 
~ 15 change, when comparing the survey count rate with 
the maximum count rate of the pointed data. 



d 




1 1.5 

CR (0.9-8.4) 



Fig. 3. Count rate CR in the soft (0.1-0.9 keV) ROSAT energy 
band versus count rate in the hard (0.9-2.4 keV) band, each 
normalized to the mean count rate in the corresponding band. 
There is correlated variability between both bands. 




Fig. 4. Count rate versus channel number in the range 0.11-2.4 
koV: Low-state data (dotted line) as compared to normalized 
high-state data (solid line). 



4.2. Spectral variability 

To check for correlated spectral changes in the hard and 
soft ROSAT energy range, we have divided the data in 
a hard and soft band, with a dividing photon energy of 
0.9 keV. This value was chosen to include most of the 
warm absorption features in the soft band, which extend 
throughout this whole region for NGC 3227 (imlike in 
other objects, where the higli-ionization oxygen and neon 
edges dominate). The soft band also includes the effects of 
the cold absorbing column, which consequently cannot be 
disentangled from that of the ionized material (see below), 
whereas the hard band is dominated by the powerlaw con- 
tinuum. We find correlated variability between both bands 
within the errors (Fig. 3), implying that no strong spectral 
changes have taken place. 

To further constrain the spectral variability during the 
observation, and to check for changes in the ionization 
state of the warm absorber in more detail, the total ob- 
servation was split (i) into individual subsets (referred to 
as 'orbits') containing about 1000 to 2000 photons each, 
time-ordered, and merged in such a way, that only pho- 
tons of comparable source flux-state are combined, (ii) 
into only two different subsets, 'high-state' (more than 
0.51 counts/s) and 'low-state'. The statistics are better in 
the latter approach, but it tends to smear out trends if 
the spectral shape at a fixed flux state depends also on 
the history. 

(i) We find constant absorption structure, described by 
the ionization parameter U, within the errors throughout 
the observation, independent of variations in the intrin- 
sic luminosity (log U = 1.05, -1.04, -1.05, -1.01, -1.01, 
-1.00, -1.02 for orbits 1 to 7, respectively and Fx = -1.9). 
The cold and warm column densities have been fixed to the 
values determined for the total observation in this analy- 
sis. The same result holds for a fiatter underlying power- 
law and the dusty models. The constancy of the absorption 
structure allows to estimate an upper limit for the den- 
sity of the warm absorber. The recombination timescale 
of the gas is given by i^ec ~ n.-^ arj.\ ; (e.g. Krolik 
& Kriss 1995), where nj/ni+i is the ion abundance ratio 
of the major coolant, cci+i^i the corresponding recombina- 
tion rate coefficient (Shull & Van Steenberg 1982), and rie 
the electron density. For the best-fit warm absorber model 
(Table 1; values in brackets refer to the dusty model) we 
find rie ^ 0.8(0.2) x 10^^t~J. No reaction of the warm 
material during the low-state at the end of the observa- 
tion implies rie ^ 7(2) x 10^ cm~^. No reaction during 
the total observation, still using tj-gc as an estimate, yields 
ric ^ 9(2.5) X 10^ cm-3. 

(ii) If U is taken to be the only free parameter, the 
same result as for the study of individual orbits is found. 
However, a more complicated sitiiation with compensating 
changes of several parameters cannot be excluded. When 
A/w is left as an additional free parameter, there is a slight 
trend for U to be lower in the high-state data (which is not 
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the way one would expect the warm absorber to behave). 
When U, Nyj and A'^h arc all left free, wo find higher U in 
the high-state data, but also larger Nu (reflecting the ten- 
dency shown in Fig. 4, i.e. a deficiency of very low energy 
photons and an excess around the location of the oxygen 
absorption edges in the high-state). However, the fitted 
quantities come with larger errors, and within these, they 
are still consistent with being constant. These trends hold 
independent of the exact value of the underlying power- 
law index, i.e. cannot be mimicked by variability of Fx 
. They are only briefly mentioned here and will not be 
discussed further. In order to study short-timescale spec- 
tral variability in detail, much better photon statistics are 
required. 
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Fig. 5. Warm-absorbed X-ray spectrum of NGC 3227 (dashed: 
dusty model) between 0.1 and 10 keV, corrected for cold ab- 
sorption. The straight line corresponds to the unabsorbed pow- 
erlaw. The horizontal line at the bottom brackets the ROSAT 
energy range. The arrow marks the position of the OVII edge. 
No iron edge is predicted to arise from the warm material. 



5. Discussion 

5.1. Evidence for a dusty warm absorber 
5.1.1. X-ray spectral fits 

Two spectral models are found to be of similar success, 
one is a dust-free warm absorber, the other one invokes 
dust mixed with the warm material. 

The presence of dust modifies the X-ray spectrum by 
absorption and scattering (e.g. Martin 1970, Martin 1978, 
Voit 1991, Laor & Draine 1993, Czerny et al. 1995). It also 
influences the conditions in the gas by several heating and 
cooling processes (e.g. Draine & Salpetcr 1979, Ferland & 
Netzer 1979, Baldwin et al. 1991, Netzer & Laor 1993), like 
gas heating by grain photoionization and grain heating by 
collisions with gas particles. The clearest signature in the 



soft X-ray region is the presence of absorption edges that 
arc produced by inner-shell photoionization (e.g. Martin 
& Rouleau 1991; their Figs 4,5). The largest edge is that 
of carbon at an energy of ~280 eV (cf. our Fig. 5); the 
shift of the edge energy due to solid state effects (Martin 
1970) is only of the order of a few eV (e.g. Greaves et al. 
1984). 

Due to the low abundance of dust as compared to the 
gas, and for typical cold Galactic columns, X-ray dust fea- 
tures are usually not detectable, since it is difficult to dis- 
tinguish from gas-phase absorption due to the same cle- 
ment if the gas is cold. However, since warm gas is more 
transparent to very soft X-rays and typically of larger col- 
umn density, dust signatures get more clearly visible. Still, 
the carbon edge itself cannot be resolved with ROSAT . Its 
presence clearly reduces the count rate in the low-energy 
channels, though. 

Since, for relatively high values of [/, dust effectively 
competes with the gas in the absorption of photons (e.g. 
Laor & Draine 1993), the spectral model involving a dusty 
absorber requires a higher ionization parameter than the 
dust-free one (Table 1). A slightly larger total column den- 
sity of the warm absorber is found for the dusty model, 
reflecting the depletion of gas-phase oxygen which is partly 
bound in the dust. 

Based on the quality of X-ray spectral fits only, no dis- 
tinction between the dust-free and dusty warm absorber 
is possible in the case of NGC 3227. Arguments in favor of 
the latter model are given below. These are derived from 
a combination of the X-ray results with optical and UV 
properties of NGC 3227. 

5.1.2. Comparison of X-ray, optical, and UV data 

As will be detailed in the following, there is evidence for a 
large amount of cold material along the line of sight to the 
active nucleus of NGC 3227 which should imprint its pres- 
ence on the soft X-ray spectrum. The evidence comes from 
emission line reddening: The broad Ha/H/? ratio as mea- 
sured by, e.g., Cohen (1983) implies an extinction of A^ ~ 
1.4 (to derive A^, we have always assumed an unreddcned 
value of Ha/H/3=3.1; if there is a very strong contribu- 
tion from collisional excitation, the estimated amount of 
reddening is reduced). Due to heating by the radiation 
field of the central continuum source, dust is not expected 
to survive in the BLR (e.g. Barvainis 1987, Netzer 1990, 
Laor & Draine 1993). Therefore, the material responsible 
for the broad line reddening should be located along the 
line of sight, between observer and BLR. A dusty outer 
BLR dominating the line emission with no clouds along 
the line of sight seems to be excluded by results from re- 
verberation mapping. These show the bulk of the emission 
to originate very near the center (r w 17 Id; Salamanca et 
al. 1994, Winge et al. 1995) where dust will be destroyed. 

If the dust that extincts the broad lines is accompa- 
nied by an amount of gas as typically found in the Galactic 
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interstellar medium, a large cold absorbing column is ex- 
pected to show up in the soft X-ray spectrum. For the 
warm absorber fit to the X-ray spectrum, we found a cold 
column of A^h = 0.55 x 10^^ cm~^, which is larger than 
the Galactic value, 7V^°' ~ 0.22 x 10^^ cm'^, but only 
corresponds to Ay=0.3. 

Mixing dust with the warm gas could supply the red- 
dening without implying a corresponding cold column. For 
the dusty warm absorber description of NGC 3227, the 
column density of the ionized gas is A^w — 10^^ ''^ cm~^ 
(Table 1), resulting in Ay ^ 3.1 which now somewhat 
overpredicts the reddening. We suggest a slightly modified 
grain size distribution with a higher abundance of larger 
dust grains, more similar to the large-R grains observed in 
the Orion nebula (e.g. Mathis & Wallenhorst 1981, Bald- 
win et al. 1991) that are characterized by less extinction 
in the optical and UV. Mixing Orion-type dust within the 
warm material, instead of the standard MRN-type dust, 
the X-ray fit remains unchanged within the error bars. It 
also has to be kept in mind that the value of Fx is not 
well constrained by the present observations. For a flatter 
powerlaw (Fx =-1.6; Ptak et al. 1994 instead of Fx =- 
1.9; George et al. 1990) the column density of the dusty 
warm material (Sect. 3.2.2) agrees with that estimated 
from broad line reddening without requiring changes in 
the dust properties. 

Other evidence for reddening along the line of sight 
is also provided by the rather steep observed UV con- 
tinuum of NGC 3227. Using the lUE measurements in 
Courvoisier & Paltani (1992; LWP 2416 and SWF 21778), 
we find a powerlaw between 2206A and 1253A of index 
aiUE — —2.9, whereas values of — —1.4 are typically 
found (e.g. Kinney et al. 1991). For comparison, extinc- 
tion with Av = 3.1 would change an intrinsic index of 
Q;iue = -1.4 to ~ —4.1 (using the extinction curve as 
listed in Osterbrock (1989; his Table 7.2)). In the con- 
text of the dusty absorber scenario, this again suggests 
dust with reduced UV-extinction properties as mentioned 
above. 

Several determinations of reddening towards NGC 
3227 are hsted in Table 2. 

The cold absorbing column of Nn = 0.55 x 10^^ cm~^ 
required for a successful spectral fit of a dusty absorber is 
still larger than the Galactic value. It is consistent with 
Mundell et al. (1995b), who find evidence for HI 21cm 
absorption towards the continuum-nucleus of NGC 3227, 
with an estimated column density of A^h — 0.55 x 10^^ 
cm~^. 

Finally, we note that even stronger reddening than that 
of the broad lines is occasionally reported for the narrow 
lines, which may, or may not, be linked to the presence of 
the warm absorber. The narrow Ha/H/3 ratio measured 
by, e.g., Cohen (1983) implies an extinction of ~ 1.2. 
Based on an observation taken at a different time (Dec. 
1993), Mundell et al. (1995a) report a much higher value. 
They find Ha/H/3=15, which corresponds to A^ ~ 4.5. 



Since under typical NLR conditions the emitted narrow 

Balmer lines shoiild always be close to their case B re- 
combination (Brocklehurst 1971) value, a change in the 
Ha/H/3 ratio has to be attributed to extinction by dust. 
However, dust intrinsic to the narrow line clouds is not 
expected to in- and decrease by factors of several within 
years, suggesting the extinction to be caused by an exter- 
nal absorber along the line of sight. The expected accom- 
panying very large cold absorbing column is not seen in 
the ROSAT X-ray spectrum. Both observations are not 
simultaneous, but they are only separated by 6^ months. 

Dust internal to the warm gas could relax the strong 
time constraint set by the quasi-simultaneous optical and 
X-ray observations; but would still require a very spe- 
cial geometry (like the NLR emission being dominated 
by clouds directly along the l.o.s. (consistent with Arribas 
& Mediavilla 1994; their Table 2), located on the far side 
only (i.e. behind the nucleus as seen from the observer), 
and absorbed by a fast-moving dusty cloud on the near 
side). The scenario is not discussed further here; a con- 
firmation of the puzzling optical observations would be 
valuable. 

5.1.3. Alternatives and implications 

An alternative possibility to explain the comparatively low 
cold column observed in the soft X-ray spectral region (as 
compared to the column inferred from broad line redden- 
ing) is that we see mainly scattered soft X-rays. In that 
case, however, the X-ray emission of NGC 3227 should not 
be that rapidly variable (Ptak et al. 1994; our Sect. 4.1). 

An alternative, that cannot be excluded, is a factor of 
a few overabundant dust in the cold column, or a strong 
dominance of small grains, to explain the observed amount 
of BLR reddening. It requires stronger deviations of the 
dust properties from Galactic MRN dust, as compared to 
the possibility of dust internal to the warm material, but 
certainly cannot be excluded with available data and we 
will include the dust- free warm absorber description of the 
data in the following discussion. 

In this context, one should re- mention the early X-ray 
observations of NGC 3227, for which the presence of large 
cold A^H was claimed. Since these had no spectral resolu- 
tion in the soft energy range, and no warm absorber model 
was applied to the data, the signature of the ionized mate- 
rial may well have mimicked a large cold absorbing column 
(note that the warm absorber modifies the intrinsic X-ray 
spectrum well above 1 keV; Fig. 5). To further test this, 
we have re-fitted the model points of the best-fit warm 
absorber model, as well as the directly observed ROSAT 
X-ray spectrum, by a powerlaw model with free cold ab- 
sorbing column. We have cut off all data points below 1 
keV. In this case, we find indeed a very large cold column 
of A^H — 3 X 10^^ cm~^, similar to what was reported for 
early X-ray observations, but which is, in fact, mimicked 
by the recovering warm-absorbed spectrum. 
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To summarize, the extinction towards NGC 3227 
seems to be caused by either (i) a warm absorber with in- 
ternal dust, or (ii) overabundant dust in the cold column. 
We tentatively favor the dusty warm absorber since it re- 
quires less strong deviations from typical Galactic dust 
properties. Simultaneous observations of emission lines, 
UV continuum, and soft X-ray properties of NGC 3227 
will be of high value. 

Table 2. Indicators of reddening/cold absorbing material to- 
wards NGC 3227. To derive from emission lines, we have 
assumed an intrinsic intensity ratio of Ha/H/3 = 3.1. The con- 
version from Av to iVn and vice versa was performed using a 
standard Galactic gas-to-dust ratio (e.g., Bohlin et al. 1978). 



method 






ref(^) 




[mag] 


[lO^Vcm^] 




broad line reddening 


1.4 


2.59 


[2] 




1.4 


2.59 


[6] 


narrow line reddening 


4.9 


9.01 


[7] 




1.2 


2.21 


[2] 




1.7 


3.13 


[3] 




4.5 


8.30 


[1] 


W/j, silicate feature 


7.5 


13.8 


[4] 


HI 21cm absorption 


0.3 


0.55 


[5] 


X-ray absorption 


0.3 


0.55 


[8] 




3.1 


5.62 


[8] 



Intrinsic to NGC 3227. Cold column density, in- 
cluding the Galactic contribution. Warm column density. 

References: [1] Mundell et al. (1995a), [2] Cohen (1983), [3] 
Gonzales Delgado & Perez (1997), [4] Lebofsky & Rieke (1979) 
and assumptions as in Reichert et al. (1985), [5] Mundell et al. 
(1995b), [6] Winge et al. (1995; their Table 4), [7] Rubin & 
Ford (1968; the authors comment on a possible underestimate 
of the H/3 flux, and estimate the uncertainty to be less than 
about 50%), [8] this work. 



5.1.4. Dusty warm absorbers in other active galaxies ? 

Notwithstanding the cautious remarks of the previous sec- 
tion, we may address the question about the possibility 
of all warm absorbers being dusty. E.g., an identification 
with the dusty transition region between BLR and NLR, 
proposed by Netzer & Laor (1993), or material evaporat- 
ing from the molecular torus, would then come to mind. 

The presence of dusty warm gas in the quasar IRAS 
13349-1-2438 was suggested by Brandt et al. (1996) 
through an analysis of the multi- wavelength properties of 
this object, and a similar model mentioned to possibly 
apply to MCG-6-30-15. On the other hand, dust internal 
to the warm gas was found to be no viable description of 



the observed X-ray spectrum of the narrow-line Seyfert 1 
galaxy NGC 4051 (Komossa & Fink 1997a). In general, we 
find the absorption structure to change towards a relative 
increase in the depths of absorption edges of more lowly 
ionized species when including dust in the warm absorber 
models (see Fig. 3 of Komossa & Fink 1997c) leading to 
an apparent Battening of the soft X-ray spectrum that 
can be partly traced back to the presence of neutral metal 
ions bound in the dust that modify the X-ray absorption 
structure (Sect. 5.1.1). One consequence is, that it is more 
difficult to reproduce apparently steep soft X-ray spectra 
(which require strong edges of highly ionized oxygen and 
neon to dominate) as, e.g., observed in narrow-line Seyfert 
Is. 

In order to assess the frequency of the presence of dust 
within warm absorbers, objects with apparently Batter 
than usual powerlaw spectra arc good candidates for fur- 
ther study, in combination with evidence for excess red- 
dening along the line of sight. The clearest signature of 
the presence of dust is the Carbon edge at 0.28 kcV. To 
uniquely detect this feature, sensitivity at this energies 
and particularly good spectral resolution is needed. The 
next mission that may achieve this is AXAF. 

5.2. Further properties of the warm absorber 

Location of the ionized material. 

(i) For the dust-free best fit (logU ~ —1-0), the 
density-scaled distance of the warm absorber is r ~ 2 x 
10^"^ /y/n7 cm, where n = lO^nr. This compares to the 
typical BLR radius of r ~ 17 Id = 4 x 10^^ cm, as de- 
termined from reverberation mapping (Salamanca et al. 
1994, Winge et al. 1995). 

(ii) In case of dust internal to the warm absorber, the 
gas density has to be less than about 10^ cm^"^, i.e. r > 8x 
IQ^^ / ^Jrvj cm, to ensure dust survival. Constant ionization 
parameter U throughout the observed low-state further 
implies n ^ 2 x 10^ cm^^ (Sect. 4.2); or even n ^ 2.5 x 
10^ cm~'^, based on the assumption of constant U during 
the total observation (less certain, due to time gaps in 
the data). The ionized material should be located at least 
between BLR and NLR; consistent with the radius of the 
BLR as determined from reverberation mapping. 

Thermal stability. The thermal stability of the warm 
material is addressed in Fig. 6. All models were calcu- 
lated with the code Cloudy. Betwcieii the low-temperature 
(T ~ 10" K) and high-temperature {T ~ 10^"^ K) branch 
of the equilibrium curve there is an intermediate region 
of midti-valucd behavior of T in dependence of U/T, i.e. 
pressure (e.g. Guilbert et al. 1983, their Fig. 1; for a gen- 
eral discussion see also Krolik et al. 1981, Netzer 1990, 
Reynolds & Fabian 1995). It is this regime, where the 
warm absorber is expected to be located, with its temper- 
ature typically found to lie around T ^ 10^ K. In case 
T is multi- valued for constant U/T, and the gradient of 



St. Komossa, H. Fink Dusty warm absorber in NGC 3227 



9 



7 



5 




-6 -5.5 -5 -4.5 

log (U/T) 

Fig. 6. Equilibrium temperature T versus U /T. The equilib- 
rium curves are shown for different ionizing continua and gas 
properties. Solid line: mean Seyfert continuum (as described 
in Sect. 3.2.1) with Fx = -1.6 (left), Fx = -1.9 (right); dot- 
ted: observed multi- wavelength continuum of NGC 4051 with 
Fx = -2.3, adopting solar abundances in all cases; dashed: Fx 
= -1.9, with depleted abundances and dust, as described in 
the text. The fat symbols mark the location of the warm ab- 
sorber in NGC 3227, for the corresponding best-fit absorber 
models. For comparison, the position of the ionized material in 
MCG-6-30-15 (arrow, solid fine; Reynolds & Fabian 1995) and 
NGC 4051 at two different epochs (crosses; Komossa & Fink 
1997a,c) is indicated. 

the equilibrium curve is positive, several phases may exist 
in pressure balance. An analysis by Reynolds & Fabian 
(1995) placed the warm absorber in MCG-6-30-15 in a 
small stable regime within this intermediate region (its 
position is marked by an arrow in Fig. 6). The ionized 
material in NGC 3227 is characterized by a comparatively 
low ionization parameter. Its location in the T versus U/T 
diagram is shown for the models that provide a success- 
ful description of the X-ray spectrum. The shape of the 
equilibrium curve is clearly altered for the model includ- 
ing dust and correspondingly depleted metal abundances, 
since both, dust and metals, affect the heating and cooling 
of the gas. Consequently, due to this strong influence of 
several parameters on the shape of the equilibrium curves 
(Fig. 6), a detailed assessment of the stability and confine- 
ment of the ionized absorber in NGC 3227, and in general, 
is still difficult in this approach. 

5.3. Predictions of warm-absorber intrinsic emission and 
absorption lines 

5.3.1. Line emission 

Since warm absorbers arc strongly matter-bounded and 
the electron temperature is high, their emissivity in emis- 
sion lines is usually expected to be less than, e.g., that 
of the BLR. However, for the high values of U typical 



for warm absorbers some optical-UV high-ionization lines 
can get very strong (e.g. Nctzer 1993, Hamann ct al. 1995, 
Komossa & Fink 1995, Reynolds & Fabian 1995, Shields et 
al. 1995), and one of the known high-ionization emission- 
line regions in active galaxies may be identified with the 
ionized absorber. E.g., in many approaches to model the 
NLR, HeIIA4686 and other high-ionization lines turned 
out to be systematically undcrprcdictcd. Shields et al. 
(1995) suggested a matter-bounded BLR component to 
act as warm absorber. 

From a sample of individual objects c;liosen to inves- 
tigate a possible warm-absorber contribution to the ob- 
served emission lines, the one in NGC 3227 represents 
the low-ionization end of the known warm absorbers. Fur- 
thermore, the optical spectrum of NGC 3227 shows broad 
wings in Ha that do not follow continuum variations (Sala- 
manca et al. 1994). Interpreting large line- width as close- 
ness to the nucleus and high density, the non- variability in- 
dicates a higher-than-usual ionized component, that may 
be identified with the; warm absorber. 

In the following we predict the emission lines expected 
to arise from the warm material, with its column den- 
sity and ionization parameter determined from the X-ray 
spectral fits. 

(i) Dust-free warm absorber: Several emission lines 
are found to be rather strong, compared to the absorber- 
intrinsic H/3 emission. The strongest arc [FeX]A6374/H/3 
= 1, HeIIA4686/H/3 — 1.1, in the optical region, and, e.g., 
CIVA1549/H/3 = 170, NVA1240/H/3 = 15, OVA1218/H/3 
= 18, OVIA1035/H/3 = 74, NeVIIIA774/H/3 = 0.3, in the 
UV - EUV. However, the absorber-intrinsic luminosity in 
H/3 itself is about a factor of 30 weaker than the mean 
de-reddened broad observed H/3 emission, Lh/s — 10^°-^ 
erg/s. Scaling the line emission of the warm absorber cor- 
respondingly, to allow a judgement on the detectability 
of these features, results in [FeX]A6374,^,a/H/3obs = 0.03, 
CIVwa/H/3obs = 6, etc. CIV is the strongest predicted line. 
Its intensity also depends on density (e.g., CIVwa/H/3obs 
= 4 for log riH = 6) and on EUV continuum shape (e.g., 
CIVwa/H/3obs = 3 for an additional EUV black body as de- 
scribed in Sect. 3.2.1). To search for a warm-absorber con- 
tribution to individual lines, high-resolution, good quality 
spectra are needed. 

(ii) Dusty absorber: In case of dust internal to the 
warm material, the overall emissivity is reduced and there 
are no strong emission lines expected to be observed in 
the IR to UV spectral region. The absorber-intrinsic ratio 
of, e.g., CIV/H/3 is weakened by a factor of more than 10, 
partially due to line destruction by grains. 

5.3.2. UV absorption 

A comparison of UV absorption lines and X-ray ab- 
sorption properties of active galaxies was performed by 
Ulrich (1988). Her analysis of an lUE spectrum of NGC 
3227 in the range 2000 - 3200 A revealed the presence 
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of MgIIA2798 absorption with an equivalent width of 
log Wx/\ ^ 2.8. The X-ray warm absorber in NGC 3227 
is of comparatively low ionization parameter. Neverthe- 
less, the column density in Mgll is very low, with log W\l\ 
= -5.7 (and considerably weaker for the dusty model or 
one with an additional EUV black body). Predictions of 
UV absorption lines that arise from the ionized material 
are made for the more highly ionized species, with col- 
umn densities of log A^c3+ = 17.5 (16.3) and log N-^i+ = 
16.8 (15.7); the numbers in brackets refer to the warm ab- 
sorber model that includes dust. (For comparison, these 
numbers change as follows, if the EUV black body compo- 
nent is again added to the incident continuum: log A^c3+ = 
17.1 (16.1) and log7VN4+ = 16.7 (15.6).) The correspond- 
ing equivalent widths are then \ogW\/\ = -2.91q;4 (~ 
S.Otg;^) for CIVA1549, and logWx/X = -3.01° j {-^^-it^^l) 
for NVA1240. The values are calculated for velocity pa- 
rameters h (Spitzer 1978) of 60 km/s, 100 km/s (for '-h'), 
and 20 km/s (for '-'). It will be worthwhile to search 
for these lines by obtaining high-quality UV spectra. For 
some examples of common UV - X-ray warm absorbers 
see Mathur et al. (e.g., 1994, 1997), Schartel et al. (1997). 

6. Summary and conclusions 

We have presented survey and pointed ROSAT PSPC ob- 
servations of the Seyfert galaxy NGC 3227. 

The temporal analysis of the data disclosed strong 
long-term variability, with a maximum change of a fac- 
tor ~ 15 in count rate (within ~ 3 years), and smaller 
changes on the timescale of days. 

The spectral analysis, possible only for the pointing, 
revealed no strong spectral variability. Two kinds of mod- 
els provide a successful description of the soft X-ray spec- 
trum, a powerlaw with soft excess, and a warm absorber 
with or without internal dust. The first is not further dis- 
cussed due to the unusually flat inferred underlying pow- 
erlaw continuum. 

Both, the dust-free and the dusty warm absorber give 
excellent fits, with values of the ionization parameter U 
and warm column density iVw of \ogU ~ —1.0, logA^^ — 
21.5, and logJ7 ~ -0.25, logA^w ^ 21.8, respectively. The 
clearest signature of the warm absorber with internal dust 
is the presence of a strong carbon edge at 0.28 keV that 
may be resolved with future X-ray missions. 

The absorber-intrinsic line emission and absorption in 
the optical UV was studied. UV absorption lines are pre- 
dicted. Significant line emission is only found for the dust- 
free description of the warm absorber, the strongest line 
being CIVA1549. However, no known emission-line region 
in NGC 3227 can be completely identified with the warm 
absorber. 

The dusty warm absorber model is shown to provide an 
explanation for the discrepancy between the low cold ab- 
sorbing column seen in the X-ray spectrum and the larger 
one inferred from emission line reddening. The density of 



the dusty warm material is constrained to at least n ^ 10^ 
cm^'^, to ensure dust survival, which translates into a dis- 
tance from the nucleus of r ^ 8 x 10^^/y^ cm and is a 
factor of several larger than the distance of the BLR as 
determined from reverberation mapping. 
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